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Abstract 

O" 

^^ , We study bb and cc production and the influence of nuclear shadowing at 

Ps| I LHC and RHIC energies. We find a significant reduction in the production 

^ I cross section of both charm and bottom at RHIC and LHC. Bound states 

O ■ such as T and J/^ are suppressed by this reduction in the charm production 

cross sections. Therefore, J/^ suppression may not be useful as a signature 

for the quark gluon plasma. 
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As lattice QCD results indicated, hadronic matter probably undergoes a phase transition 
into a deconfined quark-gluon plasma (QGP) phase, provided a temperature about 200 
MeV or an energy density ~ 2 GeV/fm^ is reached [|I|. High energy heavy- ion colliders at 
RHIC (Au-Au at 7^=200 GeV) and LHC (Pb-Pb at ^5=5.5 TeV) will provide the best 
opportunity to study the formation and evolution of quark-gluon plasma 0. As QGP exists 
only for a very short time (several fermi) in a small volume (~ 100 fm^), a direct detection 
of this state of matter is not possible. Thus various indirect signatures have to be used 
for its detection. The prominent ones among the suggested signatures are J/\l/ suppression 
1^, strangeness enhancement 0J^, dilepton and photon production [^|-|^. Although these 
signatures are proposed, it is not easy to make a concrete calculation of them for these 
experiments. Many uncertainties, both theoretical and experimental, exist which make it 
very difficult to claim the existence of the quark-gluon plasma based on the measurement of 
these signatures. At the present experiment at SPS most of the signatures can be explained 
by alternative methods without assuming the existence of quark-gluon plasma. In this paper 
we will study the impact of such non-QGP sources on the production of J/\l/ at RHIC and 
LHC energies. 

We mention here that due to the presence of nuclear, quark-gluon and hadronic medium 
a correct computation of J/\l/ suppression in ultra relativistic heavy-ion collisions (URHIC) 
becomes extremely difficult. We highlight here the major uncertainties associated with 
the calculation of J /'^ suppression. First of all one has to understand the production of 
J/\l' in nuclear collisions before addressing its suppressions. The main uncertainty lies in the 
accurate determination of the non-perturbative matrix element for J/\l/ production from the 
CC pair. Before a decade ago one used to use the colour singlet model for J/\l/ production 
at high energy hadronic collisions where it is assumed that the CC is produced directly as 
a colour singlet state which forms the J/'^ bound state [|T^. However, the colour singlet 
model failed to describe the high momenta J/'^ production data at the CDF experiment 
|TT[|. An effective QCD formalism based on non-relativistic quantum chromodynamics was 
developed where one assumes that the CC is formed as a colour octet state which finally 



emits a soft gluon to form a colour singlet J/\l/ |T^. This model (known as colour octet 



model) is very successful in explaining the J/\& production data at the CDF experiment at 



Fermi Lab [|TI|. In this model the non-perturbative matrix element for J/\E' formation from 
CC is determined from the experiment. However, it is not at all clear how to obtain this non- 
perturbative matrix element for high energy nucleus-nucleus collisions. The complication 
arises due to the presence of a nuclear/quark-gluon medium. The second uncertainty lies 
in the accurate determination of the nuclear partonic structure function used to compute 
CC . Once these two issues are solved then one can compute J/\l/ production in high energy 
heavy-ion collisions. We will address here how significantly the heavy quark production 
cross section is changed when modified nuclear partonic structure functions are taken into 
account. This will have a direct impact on the quarkonium production cross section in high 
energy heavy-ion collisions at RHIC and LHC. 

Before proceeding to our calculation we briefiy mention some of the other effects which 
must be considered for a correct calculation of J/\l/ suppression. First of all we mention 
that the space-time evolution of quark-gluon plasma in ultra relativistic heavy-ion collisions 
is described in various stages, namely; 1) Pre-equilibrium, 2) Equilibrium and 3) Hadroniza- 
tion. Once CC is formed, it travels in pre-equilibrium stage of the quark-gluon plasma. The 



charm quarks scatter with the hght quarks and gluons in the medium and hence change the 
energy, momentum and separation of the CC pairs. It is not easy to keep track of all the in- 
teractions a charm quark encounters during its propagation in the non-equilibrium medium. 
This is because a correct description of non-equilibrium space-time evolution of the quarks 
and gluons formed in nuclear collisions is not solved accurately. On the other hand, if a J/^ 
is formed, then its interaction with quarks and gluons in non-equilibrium medium has also 
to be taken into account. In the equilibrium stage one can argue that they are screened. But 
a real dynamical screening argument has to be given. The screening argument given in 
is valid for a static plasma as calculations were done in lattice QCD which assumes a static 
temperature. In heavy-ion collisions the temperature is local and space-time dependent and 
a screening calculation involving local space-time dependent temperature is not available. 
The J/\l/ also interacts with the hadrons in the hadronic phase |1^,|T^. Hence all these 



effects have to be considered in the study of J/\l/ suppression at RHIC and LHC. 

In this paper we address one of the above issues, i.e. we compute heavy-quark production 

at RHIC and LHC by taking modifications of the nuclear partonic distribution functions 

into account. The modification of nuclear parton distribution functions {i.e, F^{x,Q'^) ^ 

A F^i^XjQ"^) and G^(x,Q^) 7^ A G'^{x,Q'^)) has been observed in various experiments 

T5| . Here F2 (x, Q^) and G^^^\x,Q'^) are the quark and gluon structure functions of 



the nucleus (nucleon), x is the momentum fraction of the parton, and Q^ is the momentum 
transfer in deep inelastic lepton nucleus (nucleon) scattering. We take the modified parton 
distribution functions from a previous study [jl6|, where it has been shown that there is a 
strong shadowing of the parton distribution functions at low x. The idea of parton fusion 
was first formulated in [l^ and later proven in |jl8| to appear when the total transverse 



size 1/Q of the partons in a nucleon becomes larger than the proton radius to yield a 
transverse overlapping within a unit of rapidity xG{x) > Q^R^. The fusion correction in 
the free nucleon turns out to be significant only for unusually small values of x and Q^. 
However, in the case of heavy nuclei of A ~ 200, the parton recombination is strongly 
increased (see |]TB|,irU|,ED[). This strong shadowing is expected to reduce the heavy quark 



production significantly. We incorporate this feature of gluon shadowing in the calculation 
of CC production at RHIC and LHC heavy-ion collisions. 

We use perturbative QCD methods to compute charm and bottom production at RHIC 
and LHC nuclear collisions. As we are interested in comparing the results of heavy-quark 
production using unshadowed and shadowed parton distributions we will not do a next- 
to-leading order calculation for QQ production. For our purpose we consider the leading 
order gluon fusion process g{pi) + g{p2) -^ QiPs) + Q{Pi) for heavy quark-antiquark pair 
production. For simplicity a K factor ~ 2 can be multiplied with this leading order cross 
section to consider higher order processes [ETI] . The gluon fusion process is dominant at high 



energy hadron-hadron collisions ||2^. The squared amplitude for the above leading order 
process is given by 
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where M is the mass of the heavy quark Q. The Mandelstam variables s, t, u are defined by 

s = {pi+ P2)^ i = {pi - pa)^ M = (pi - p^y 

with the four- momenta pi, p2 and ps, pi of the incident and the scattered particles respec- 
tively. Using the amplitude for the 2^2 process, the partonic cross section can be written 
as 



da 1 rM^-s{l-^l-AMys)/2 

ais) = / dt— = -— ^ / dtY,\M\\ 2 

at 167rs^ 7Af2-s(i+^i-4Af2/s)/2 

Assuming Q^ = s ||23|, the a((5^) can be taken out of the t integration to find 
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This cross section in the partonic level can be used to compute the QQ production cross 
section in hadronic collisions by using the parton distribution functions inside the hadron. 
The total cross section for QQ production per nucleon in AA collision is: 



(^(s)qq= ,dxi dx2gf{xi,Q^)gf{x2,Q^)agg_,QQ{s), (4) 

J4M^/s J4M-' /(sx-i) 
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/4M2/S JAW^Ksxi) 

where s = X1X2S, and g^{x, Q^) is the gluon structure function of a bound nucleon inside 
the nucleus A. On these equations nuclear shadowing effects in g^{x, Q^) will have an influ- 
ence, as the nuclear structure functions are reduced in comparison to free nucleon structure 
functions at low x. We use the shadowing version of parton distribution functions in a 
bound nucleon from ref. |[T^ (with A = 208) together with the GRV98 set of gluon distribu- 
tion functions for a free nucleon ||2^. The strong coupling constant a{s) used in the above 
equation is also taken from |21|. 



We mention here that our calculation is based on primary hard collisions of partons. 
These primary hard collisions also produce a huge number of jets and minijets at RHIC 
and LHC energies. Unlike in pp collisions these minijets suffer from secondary collisions 
because they are large in number. Secondary collisions among these primarily produced 
minijets bring the system towards equilibrium. The secondary collisions among these mini- 
jets can also produce additional CC pairs. However, it should be checked whether this 
heavy quark production is small or large due to their large masses. It can be argued that 
the secondary charm production is very small for an equilibrium quark-gluon plasma be- 
cause the production rate behaves like ~ e~*^/"^. However, in the early pre-equilibrium 
stage of the quark-gluon plasma, the average energy carried by a parton is large and hence 
their secondary collisions can produce signiflcant secondary CC pairs. One needs a real 
calculation by using non-equilibrium distribution of the minijet plasma at RHIC and LHC 
to see whether the secondary CC pair production becomes larger/smaller than the primary 



charm production. The secondary charm production calculation is very important because 
a huge production of secondary charm may suggest an enhancement of J/\& at RHIC and 
LHC rather than its suppression. In any case, the shadowing will also decrease the minijet 
production |T^. In such situations the minijet temperature will also be reduced which in 



turn will decrease secondary charm production. Hence as far as shadowing of initial parton 
distributions is concerned, it decreases both primary and secondary charm productions. In 
the following we will present results of the cross section for heavy quark-antiquark pair pro- 
duction from primary partonic collisions at different center of mass energies to observe the 
effect of shadowing on quarkonium production. 

In Fig.l we present the cross section (in mb) for heavy quark/anti-quark pair production 
at RHIC energy as a function of the heavy-quark mass. The dotted line is obtained by using 
no-shadowing gluon distribution functions, i.e, we have used g^i^x^Q"^) = g'^{x,Q'^), with 
g'^{x, Q"^) taken from the GRV98 parametrization pi) (Note that g^{x, Q"^) used in Eq. (H) 
is the distribution function of a bound nucleon inside a nucleus A). It can be mentioned 
that our results are with respect to the nucleon-nucleon scattering in nuclear collisions. 
To estimate the total cross section in nucleus-nucleus collisions, one has to multiply our 
results by A^. The solid line is obtained by using strong shadowing of gluons at high energy 
| |16| |. It can be seen that the shadowing in gluon distribution function decreases the heavy 
quark production cross section significantly. The reduction of charm quark production cross 
section at RHIC and LHC is interesting because the suppression of J/\l/ may be due to the 
shadowing of the initial gluon distribution functions inside the nucleus rather than existence 
of QGP. However, before making any conclusions, one has to compute J/\l/ production by 
using a suitable model for hadronization. In addition to this, one has to compute secondary 
charm quarks which are produced from the rescattering of minijets being formed in high 
energy heavy-ion collisions. We have plotted the ratio of heavy-quark production cross 
sections with and without shadowing in the same figure. This is represented by the dashed 
line. Our results cover both charm and bottom quark production at RHIC. In future RHIC 
will provide an opportunity to measure open charm production. 

In Fig. 2 we present the cross section for heavy quark/anti-quark pair production at LHC 
energy as a function of the heavy-quark mass. We present our results up to Mq ~ 200 GeV 
which covers charm, bottom and top quark production at LHC. Dotted, dashed and solid 
lines correspond to the same situation as in Fig.l. The difference in charm and bottom 
quark production at LHC is found to be much larger than that at RHIC. It can be seen 
from the figure that the shadowing effect on charm production is much larger than that on 
bottom production. The ratio is close to 1 for the production of the top quark. 

In Fig. 3 we present the CC production cross section as a function of i/i. The dotted 
line is for the case without shadowing, and the solid line is with shadowing of gluons. It can 
be seen that the CC production cross section at RHIC for both cases is not much different. 
However, the difference becomes quite large when the calculation is repeated for very high 
energy nuclear collisions. At LHC energy the reduction in the cross section is quite large, 
about a factor of 10. The same is true for BB production at different energies, see Fig.4. 
We mention here that our results are obtained from the leading order partonic scattering 
processes. Our results can be multiplied by a K factor (~ 2) to take higher order processes 
into account. 

To summarize, we have computed heavy quark/ anti-quark pair production cross sections 



per nucleon-nucleon collision in RHIC and LHC heavy-ion collisions with shadowing of 
parton structure functions taken into account. It is shown that the shadowing of parton 
distribution functions significantly decreases the production cross section of heavy quarks at 
LHC. At RHIC energy this difference is not very large but still significant. One can obtain 
the total charm production at RHIC and LHC from these cross sections via the relation 
A^ = T(0)o", where the Glauber geometrical factor T(0) ~ ^;r^ for the impact parameter 
b=0 fm. This geometrical factor represents the total number of binary nucleon-nucleon 
collisions in central AA collisions, which behaves like A^^^ p5[. Although the difference of 



the bb and cc cross sections for shadowing and no shadowing at RHIC is not large, there is a 
significant difference in the total number of charm and bottom quarks for the two versions of 
shadowing. This has a direct impact on quarkonium and open charm production at RHIC. 
We suppose that a reduction of the charm quark cross section due to shadowing is very 
significant because J/\l/ suppression might occur due to the depletion of the initial gluon 
distribution functions rather than due to the existence of a quark-gluon plasma or due to 
the presence of hadronic medium. 
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FIGURE CAPTIONS 

Fig. 1. The cross section for lieavy-quark production via the parton fusion process 
99 ~^ QQ ^-t RHIC. The sohd (dotted) hne represents the heavy-quark production cross 
section (in inh) with (without) gluon shadowing (explained in the text). The dashed hne 
represents the ratio of the cross sections with and without shadowing. 

Fig. 2. The cross section for the heavy-quark production via the parton fusion process 
99 ~^ QQ at LHC. The sohd (dotted) hne represents the heavy-quark production cross 
section (in mh) with (without) gluon shadowing (explained in the text). The dashed line 
represents the ratio of the cross sections with and without shadowing. It can be seen that 
there is a strong shadowing effect on charm production (about a factor of 10, (see the ratio 
curve and discussion in the text)). 

Fig. 3. The cross section in mh from the process gg -^ CC. Solid (dotted) line is with 
(without) gluon shadowing (explained in the text). Mc = 1.5 GeV is used in the calculation. 

Fig. 4. The cross section in mh from the process gg -^ BB. Solid (dotted) line is with 
(without) gluon shadowing (explained in the text). Mb = 5.4 GeV is used in the calculation. 
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